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ABSTRACT 

Fourier power spectra and phases of a signal from a large number of radiating sources 
orbiting around a black hole are investigated. It is assumed that the individual sources 
(bright spots) are located in an accretion disc and their lifetime exceeds the corre¬ 
sponding orbital period. This model is relevant for the short-time X-ray variability 
of active galactic nuclei. Previous works on this subject were mostly concentrated on 
temporal characteristics and power spectra of observed light curves. In our present 
contribution, Fourier phases are brought into consideration and studied systemati¬ 
cally for a broad range of input parameters. In particular, conditions for the phase 
coherence are discussed. It is shown that one can discriminate between the two classes 
of models which are currently under consideration—orbital motion of a large number 
of sources versus short-lived independent flares—although parameters of the model 
are not completely arbitrary. It is also shown that predicted power spectra depend 
rather strongly on the spot distribution across the disk surface. We conclude that the 
orbital motion of the spots cannot be the only reason for the source fluctuations, but 
it certainly influences observational properties of the source intrinsic variability. 

Key words: galaxies: active - accretion, accretion discs - X-rays: galaxies - black 
hole physics 


1 INTRODUCTION 

Observational effort of recent decades brought an extensive 
amount of data abont the variability properties of active 
galactic nnclei (AGNs) in the range of wavelength from ra¬ 
dio to hard X-rays (Duschl, Wagner & Camenzind 1991; 
Miller & Wiita 1991). It is widely presnmed, although not 
yet proven, that massive black holes with accretion discs re¬ 
side in cores of AGNs (Blandford, Netzer & Woltjer 1990). 
Within this model which is also accepted in the present con¬ 
tribution, most of the X-rays originate inside inner regions 
of AGNs, a few or a few tens of gravitational radii from the 
center. The observed signal shows irregular, featureless fluc¬ 
tuations with a very complex behaviour at frequency about 
10“®-10“^ Hz. The signal has traditionally been analyzed 
by statistical approaches, the Fourier analysis in particular 
(Feigelson & Babu 1992). In the above-mentioned restricted 
range of frequency, power-spectrum density of the fluctuat¬ 
ing signal can be represented by a power-law: P{uj) oc u}~°‘ 
with 1 < a < 2 (Lawrence et al. 1987; McHardy & Cz¬ 
erny 1987; Mushotzky, Done & Pounds 1993; Lawrence & 

* Also at Scuola Internazionale Superiore di Studi Avanzati, Tri¬ 
este; Department of Astronomy and Astrophysics, Goteborg Uni¬ 
versity and Chalmers University of Technology, Goteborg 


Papadakis 1993). This form of the power spectrum can be 
accommodated by several physical models of the central re¬ 
gion which have been proposed in the literature, restricting 
partially the parameter space of each of them. However, the 
power-spectrum itself does not allow us to discriminate be¬ 
tween intrinsically different models or exclude some of the 
proposed scenarios from further consideration. Models of 
the AGN X-ray variability refer to various instabilities in 
the inner parts of the accretion disc, the disc corona or the 
jet which probably operate in a nonlinear regime (Krolik 
et al. 1991; Mineshige, Ouchi & Nishimori 1994; Ipser 1994; 
Kanetake, Takeuti & Fukue 1995). Another viable approach 
considers the observed variability as a direct consequence of 
the orbital motion of numerous irregularities (regions of en¬ 
hanced emissivity, or bright spots) on the surface of the ac¬ 
cretion disc (Abramowicz, Bao, Lanza & Zhang 1991; Wiita, 
Miller, Carini & Rosen 1991). The latter scenario has been 
formulated on a phenomenological level. The origin of the 
spots remains an open question and the mechanism, by its 
definition, turns out to be relevant only if the lifetime of in¬ 
dividual irregularities is longer than the local orbital period. 
In this respect, the possible existence of vortices (Abramow¬ 
icz, Lanza, Spiegel & Szuszkiewicz 1992) and spiral waves 
(Chakrabarti & Wiita 1994) in accretion discs has been dis¬ 
cussed in the literature. The bright-spot model has been also 
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applied to to the UV microvariability of AGNs (Mangalam 
& Wiita 1993). 

In the present contribution, we wish to clarify whether 
there is enough useful information in the Fourier spectra 
to indicate that the observed fluctuations are due to or¬ 
bital motion of separate (mutually unrelated) sources of 
radiation rather than virtually random flares. It has been 
stressed repeatedly that an appropriate approach to data 
analysis should take the physical nature of the system into 
account (cf. Vio et al. for references). It appears conceiv¬ 
able, however, that one can discriminate between the two 
alternatives—orbital motion of sources each of which sur¬ 
vives many orbital periods versus short-lived flares—even 
when parameters of the model are not fully determined. The 
reason for this hope is apparent: orbital motion of the source 
modifies the intensity of observed radiation periodically. The 
relative change of the radiation flux (count rate) is particu¬ 
larly large when the disc is seen edge-on. This effect is mainly 
due to gravitational focusing (Bao 1992; Karas, Vokrouh- 
licky & Polnarev 1992, and references cited therein), and 
it turns out that light rays must be calculated properly by 
solving the geodesic equation when inclination of the ob¬ 
server is greater than about 70°; otherwise, simulated light- 
curve profiles are wrong. On the other hand, the hope of 
finding a signature of the orbital motion in fluctuating light 
curves may be false when the sources are spread through a 
large range of radii from the black hole, so that their orbital 
periods differ too much. Also, the harmonic content of the 
light-curve profiles is then large. These facts are well known 
from studies of the power-spectra (McHardy 1989; Zhang & 
Bao 1991) and indicate that the bright-spot model should be 
considered as one of viable explanations of the quasi-periodic 
oscillations (see, e.g., Miller & Park [1995] for a recent list 
of references). 

In the next section, properties of the Fourier transform 
relevant for the present work are briefly summarized. Fourier 
phases are discussed because it turns out that power spectra 
do not contain enough information. Details of our model of 
orbiting sources are then specified and results of a systematic 
exploration of its parameter space are illustrated. 


2 FOURIER POWER SPECTRA AND PHASES 

Power-spectrum analysis is useful for detecting periodicities 
present in a signal, and for this reason it has been discussed 
extensively in the astronomical literature and applied to var¬ 
ious problems. The power spectrum, however, contains less 
information than original data or their Fourier transform 
from which the power spectrum has been constructed. The 
missing portion of information can be reconstructed from the 
Fourier phase. Physical interpretation of the Fourier phase 
and its relation to properties of the corresponding system 
is much less understood than the power spectrum. Analy¬ 
sis of simple model cases suggests that it may be possible to 
reveal hidden periodicities by studying Fourier phases. Lanc- 
zos & Gellai (1975) performed the Fourier analysis of ran¬ 
dom, computer-generated sequences and noticed patterns in 
the distribution of the Fourier phases. Their aim was to in¬ 
vestigate periodicities which are present in pseudo-random 
numbers. (A historical note which is added to the report 
by Lanczos & Gellai suggests that the original motivation 


Figure 1. Fourier analysis of the signal from a single orbiting 
spot is presented as an illustrative example which helps to under¬ 
stand more complicated situations with a large number of spots. 
No noise component was added to the signal. Parameters of this 
figure are summarized in the first row of Tab. ^ The counting rate 
(a) is mean subtracted and normalized. Notice how the phase co¬ 
herence of the signal is reflected in panels (d)-(f). 

to consider Fourier analysis of random sequences stemmed 
from Einstein’s [1915] article.) Non-uniform distribution of 
phases and their clustering to certain angles is usually re¬ 
ferred to as a phase coherence. Krolik, Done & Madejski 
(1993) have explored the phase coherence of X-ray light 
curves of five Seyfert galaxies observed by EXOSAT but 
found no statistically significant coherence in their data. (It 
has been demonstrated later that the EXOSAT data for one 
of their sources, then a popular candidate for the AGN pe¬ 
riodic variability, were dominated by a cataclysmic variable 
in the Galaxy; cf. Madejski, Done, Turner, Mushotzky et al. 
1994.) 

2.1 Details of the model 

In this paragraph we specify the bright-spot model with suf¬ 
ficient generality and ask whether its relevance for under¬ 
standing AGN variability can be tested by supplementing 
discussion of power spectra with analysis of Fourier phases. 
The mass of the central black hole in AGNs is usually es¬ 
timated to be in the range of M « 10®-10^^ Mq. The 
mass of the accretion disc is neglected and gravitational 
field is described by a vacuum spacetime of a nonrotat¬ 
ing (Schwarzschild) black hole (Misner, Thorne & Wheeler 
1973). The corresponding gravitational radius of a black hole 
is Rg = 2c~^GM « 10“®M8 pc, and characteristic time 
tg = 2c~^GM ~ 10®Ms sec, where Ms = M Mq). Ge- 
ometrized units, c = G = 1, are used hereafter. All lengths 
and times are made dimensionless by expressing them in 
units of M. 

It is assumed that the disc is geometrically thin and ro- 
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tates with the Keplerian circular frequency. A generalization 
to geometrically thick discs is straightforward once the disc 
height and its angular velocity are specihed as functions of 
radius. The main difference between geometrically thick and 
geometrically thin discs consists in a possibility of eclipses 
when observer inclination is large, but we neglect these ad¬ 
ditional free parameters of the model in the present contri¬ 
bution; it turns out that the restriction is irrelevant to the 
question whether resulting light curves are phase coherent or 
not. Profiles of frequency-integrated light curves of an orbit¬ 
ing source of light were originally studied by Cunningham & 
Bardeen (1973). Later, a number of authors employed var¬ 
ious computational schemes which are needed for efficient 
evaluation of observed radiation flux resulting as a super¬ 
position of contributions from a large number of individual 
sources. In the present work, an analytical htting-formula 
for the observed counting rate Ra, Oq) from each of the 
spots (in arbitrary units) was employed: 


F{^,Ra,eo) = - 1 )'/") 


cos 9o 


1 -I- sin [27r(i/j -|- P 4 - |) + f ] 

-I- cos“^ Bo (pi -I- peR^J^) exp [—P 222 ], (1) 


Zl = P 5 COS^^^ 60 -\- ps COS®^^ 60 , 

_ I 1 |9/5 
Z2 = \p-^\ ; 


Rs is the radius of the orbit of the source. Bo is observer 
inclination (0o = 90° corresponds to the edge-on view of 
the disk), and ip oc is the normalized phase of the 

periodic signal (0 < <p < 1). Eq. (1) takes into account 
the fact that observed energy of radiation is different from 
radiation energy in the local rest frame of the source of light 
(Doppler effect, gravitational redshift). Focusing of light by 
the central black hole is also included. 

Several restrictions have been imposed in derivation of 
the fitting-formula (1) (Karas 1996): (i) The bright spot 
is located in a Keplerian circular orbit in the equatorial 
plane of a black hole; (ii) Radius of the orbit satishes condi¬ 
tion 3Rg < Rs <45 Rg (iii) Local emissivity of the bright 
spot is isotropic (in the rest frame comoving with the spot) 
and it decreases exponentially from the center of the spot; 
(iv) Characteristic size d of the spot (defined as a distance 
from the spot center to the point where emissivity decreases 
by factor 1/e with respect to its value in the spot center) 
satisfies d <C 7?^; (v) Inclination angle satishes condition 
20° < Bo < 80°. 

The numerical code which produced simulated light 
curves was described by Karas, Vokrouhlicky & Polnarev 
(1992), while derivation and further details concerning 
eq. (1) have been discussed in the above-mentioned paper 
(Karas 1996). Here, we only give the values of parame¬ 
ters Pk which were obtained by non-linear least-square ht- 
ting to numerically simulated light curves: pi = 0.021696, 
P2 = 190.7236, P3 = 0.3476, p4 = -0.0018, ps = 3.5106, 
P6 = —3.6 X 10“®, P 7 = 0.0124, Ps = —0.0231. Only a non¬ 
rotating black hole is considered in the present contribution 
because parameters Pk are not sensitive to the black-hole an¬ 
gular momentum. It is worth noting that the analytical form 
of eq. (1) proves to be more ehicient than previous numerical 


Figure 2. As in Fig. ^ but with 100 spots orbiting in the range 
of radii 6 Rg < Rs < 2b Rg. Low values of a and p (Tab. indi¬ 
cate the phase-coherence which can also be inferred from varying 
histograms in the panel (d), although the effect is less apparent 
than in Fig. |^. 


Figure 3. As in Fig. ^but with 5 % noise component added. One 
observes a decrease in the phase coherence due to the noise. See 
the text for a detailed description. 


approaches, while preserving the characteristic shape of the 
light curve to which gravitational focusing and the Doppler 
effect both contribute. We found no indication that our re¬ 
sults are affected by the above-given assumptions (i)-(v) 
which were needed only in derivation of the htting-formula 
(1) from simulated data by the least-square procedure. 

Individual contributions to the total hux were gener¬ 
ated with random orbital phases and the resulting flux was 
submitted to the Fourier analysis. The light-curve prohle 
from a single source (and without a noise) is, naturally, pe- 
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riodic in time. The idea of these models is that a superpo¬ 
sition of many contributions with different periods results 
in a fluctuating signal. According to the present status of 
the bright-spot model, it appears that the fluctuations due 
to the orbital motion cannot be the only reason of AGNs 
variability but they are certainly coupled with intrinsic fluc¬ 
tuations and influence the resulting observational properties 
of AGN. Denoting the total flux from Ns sources by 

Ns 

F{t,eo) = '^F„{t,Rs,9o), (2) 

n=l 


and carrying out the Fourier transform, F{t) —> F{u!), one 
obtains for the Fourier power spectrum and the Fourier 
phase (Born & Wolf 1964): 


Pico) 


1 

2 


lim 

t —>oo 


\f{uj) 


(3) 


())(cu) = arctan 


»e[F(tj)] 


(4) 


The phase (f) takes values from interval (—7r,7r). The two 
boundary values, ±7r, are identified with each other, and 
play thus no particular role. It is then possible to smooth the 
power spectrum and, within a restricted range of frequency, 
to approximate the resulting curve by a power-law, 


P{uj)=uJ-°‘. 


(5) 


The phase coherence is expected between harmonics cor¬ 
responding to the spots located at the same radius, but 
it is difficult to estimate how the coherence will be sup¬ 
pressed when a range of radii is considered. A large number 
of model parameters can conveniently be investigated by 
checking whether a is within the expected interval of values 
and then, for selected cases, performing a formal test for the 
phase coherence. 

We define an average relative distance of M phase 
points by 

M 

(6) 

m=2 

with 


= Min{|<))^ - - ferrl} , fc =-1,0,1. 

Frequency is assumed to fall within a limited interval 
{^min = l/Tmax, ^max = l/Tmin) where Prnin = ^(T^niin) 
and Tmax = r(i?max) are, respectively, orbital periods at 
the minimum and the maximum radius of the spot distribu¬ 
tion. M, the total number of the phase points in the spec¬ 
trum, can, in principle, acquire an arbitrary large value in 
our tests, but it is very restricted in real observations be¬ 
cause of their temporal sampling. Apparently, 0 < cr < 1, 
with values near zero indicating the kind of coherence where 
the phase changes only slowly while values near unity indi¬ 
cate complete incoherence of the phase points. 

Another, also rather crude test of the phase coherence is 
taken from §3.2 of Krolik et al. (1993). The frequency-phase 
plane is binned and the consistency of the distribution with a 
uniform distribution is tested by calculating statistic and 
corresponding probability p. The binning must be chosen 
appropriately so that there are enough phase points in each 


of the bins, Mi. We accepted a uniform binning into two 
frequency bins between cUmin and cUmax, and four bins in 
phase angles. Thus we write, in standard notation, 

= p = r(|,ix"). (7) 

i = l 

Let us note that the two tests described above are not equiv¬ 
alent in general, but they are both capable of indicating the 
phase coherence which occurs, for example, in a simple sine- 
type signal. Since no formal test of the phase coherence with 
a general validity is available, a direct inspection of graphs 
showing the phases also helps. 

We studied a simplified model in which individual 
sources are identified with bright spots. They are distributed 
across the disc surface with the number density and maxi¬ 
mum observed fluxes being given by power-laws: 

n{R) oc oc R-°‘*. (8) 

In addition, we assumed that the spots have a fixed radius 
and are destroyed after several revolutions. We checked that 
results are not very sensitive to the latter two assumptions 
about the size and the lifetime of the spots; it is only re¬ 
quired that the size d satisfies relation 0 < d Rs and 
the lifetime ta 3> Rs ■ However, power-spectrum indexes 
are rather sensitive to the spots distribution, as discussed 
later. Our work extends previous calculations in which in¬ 
formation in Fourier phases was not taken into account and 
the light-curve profiles of individual spots were either ap¬ 
proximated by a simple box-type pulse (Zhang & Bao 1991; 
Abramowicz et al. 1991) or their parameters were assumed 
to fall in a more restricted range (Bao 1992). 


2.2 Results 

Description of the spots’ characteristics in terms of power- 
laws represents a simplified model. We carried out a system¬ 
atic investigation of the parameter space of this model. Our 
results are illustrated in Figures First, it is useful to 
consider a single spot because in this case the correspond¬ 
ing light curve has a j^rticularly simple form. The situation 
is explored in Figure hi which consists of six panels: 

(a) Light curve in the time domain. The curve is mean 
subtracted and normalized to the maximum counting rate. 
A detail of the curve which covers several periods is shown. 
Naturally, the light curve of a single spot is strictly periodic. 

(b) Power spectrum P{ui) in the log-log plane. Frequency 
range is chosen to cover (cUmin, Wmax) (a small overlap is nec¬ 
essary because, with only one spot, boundary frequencies de¬ 
generate to Wmin = Wmax and Rs = Rmin = Rmax). The main 
peak in the power spectrum corresponds to the fundamental 
frequency while the smaller peaks are harmonics. 

(c) Fourier transform in the complex plane. While this 
graph often exhibits interesting patterns, it is difficult to in¬ 
terpret. Gircles correspond to frequency points to < Wmax; 
crosses correspond to higher frequencies (they are concen¬ 
trated with a small amplitude near the origin). 

(d) Histogram of the Fourier phases, (^(w). The graph 
is normalized to unity and shows a relative distribution 
which can be compared with other graphs representing dif¬ 
ferent numbers of the phase points. Unequal distribution 
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Table 1. Fourier analysis of the signal generated by bright spots orbiting around a Schwarzschild black hole. Selected cases are illustrated 
in corresponding Figures, as indicated in the first column of the Table. 


No noise 5 % noise component added 


Fig. 

Ns 

-^min 

-^max 

do 


ai 

a 

a 


P 

a 

(T 

X^ 

P 

Q 



1 

56 

56 

25 

- 

- 

- 

0.021 

107.8 

0 

- 

0.388 

30.55 

0.0050 

- 


5 

100 

6 

25 

50 

0 

1 

3.0 

0.485 

662.2 

0 

2.9 

0.607 

210.96 

0 

0.3 



200 

6 

25 

20 

0 

1 

3.2 

0.537 

5.77 

0.567 

3.1 

0.640 

5.01 

0.659 

0.3 


- 

200 

10 

30 

80 

0 

0 

1.8 

0.611 

9.87 

0.196 

1.8 

0.619 

2.84 

0.899 

1.7 

- 

200 

6 

25 

80 

1 

0 

1.7 

0.600 

52.89 

4 X 10 -® 

1.7 

0.640 

20.03 

0.0055 

1.0 

- 

200 

30 

60 

80 

-1 

0 

2.1 

0.615 

9.84 

0.198 

2.1 

0.654 

15.46 

0.0305 

2.3 

- 

200 

6 

25 

80 

-1 

0 

1.7 

0.600 

52.89 

4 X 10 -® 

1.7 

0.640 

20.03 

0.0055 

2.3 


Ns is the number of spots in the model; (iJmini ^max) is the range of radii; Bq is observer inclination: an, a^ are parameters specifying 
distribution of the spots (eq. |^); a and a are the power-spectrum indexes (eqs. and 1^); cr (eq. ^), s-iid P (stf- @1) characterize 

the phase properties of the model. 


indicates that the phase points form clusters around certain 
frequencies. On the other hand, phase incoherence yields a 
nniform distribution of the histograms. Frequency interval 
(oimin, Wmax) IS again Considered. Recall that points (;i± 180° 
are to be identified. 

(e) Cumulative distribution of the Fourier phases. This 
graph contains the same information as the graph (d), but it 
can be directly compared with Fig. 4 of Krolik et al. (1993). 
A signihcant deviation of the cumulative distribution from 
the diagonal line (dashed) indicates the phase coherence and 
can be formally tested by the Kolmogorov-Smirnov statistic. 

(f) Graph of individual phase points as a function of fre¬ 

quency. Phase coherence is visible in this graph where the 
points form several evident sequences. _ 


Fi^. y shows a pure light curve (nu additional noise). 
In order to simulate a high-frequency noise, either inherent 
to the signal or due to instrumental errors, we also super¬ 
posed simulated data with the noise (with zero mean and a 
small, typically 5% standard deviation). Typical situations 
are illustrated in Figures (where the phase coherence is 
present) and Fig. ^ (where no phase coherence is indicated). 
Parameters of these figures are given in Table □. 

The first row of Tab. corresponds to Fig.^ while sub¬ 
sequent rows describe situations with much higher numbers 
of spots which, together with a noisy contribution, form 
the predicted fluctuating signal of the source. The power- 
spectrum index, a (eq. |5|), is also calculated. The value of 
Q (the slope of the power spectrum as predicted by the cur¬ 
rent work) can be compared with a corresponding quantity 
a in the approximation of Zhang & Bao (1991): 


d = I (4 - 2Qi - On) - 1, (9) 

which is obtained from their eq. (3.14) by setting an = 3/2 
and at = 0. Eq. (^) for a has been derived analytically 
but with a number of simplifying assumptions about the 
light-curve profiles. On the other hand, corresponding val¬ 
ues of the power-spectrum slope a follow from our numer¬ 
ically derived fitting-formula (1). Let us note that a does 
not depend on observer inclination while our a shows this 
dependence. The mean value of the slope is typically in the 
range 1 < {a) < 4, but {a) < 2 whenever the inclination 
is restricted to 0o ^ 75°. Power spectra become flatter at 
higher frequencies. Power-spectrum index is rather sensitive 
to the spot distribution as characterized by (q„ = 0 
for the uniform distribution). A general trend of a being 


anti-correlated with q„ (Eq. y) remains preserved in our 
analysis. One needs to be careful in making direct compar¬ 
isons between a and a because we always considered very 
restricted range of Rs < Rmax (presumably a reasonable 
restriction for studying the X-ray variability). In addition 
there are further, more subtle differences betwen previous 
works and our present approach: the fitting formula (1), fi¬ 
nite sizes of the spots, and a numerical package we have 
employed (Matlab v. 4.0 with the usage of its built-in pro¬ 
cedures) but we do not expect these differences to affect 
the results significantly. (Indeed, Tab. 1 only briefly sum¬ 
marizes typical results of the systematic search over the pa¬ 
rameter space; lengthy tabular material can be found at 
our Web-site, “|ittp://otokar.troja.mfi’.cuni.cz/user/karas/ 


au_www/karas/papers.htm'’.) We currently work on a de¬ 


tailed comparative analysis of different approaches which 
should determine the range of applicability of the simplified 
analytic formula (H). 


3 CONCLUSION 

Within the considered range of frequency, Hz, ob¬ 

servations indicate 1 < a < 2 and, very preliminary, no 
phase coherence in the Fourier spectra of the signal from 
AGNs. We have seen that the bright-spot model is restricted 
(but not excluded) by these facts. It is suggested that the 
procedure described above can be used to contrast theoret¬ 
ical models of AGNs variability with observations. These 
models are testable. Our present work was motivated by the 
bright-spot model of X-ray featureless variability in its orig¬ 
inal formulation with a random distribution of individual 
spots. One might speculate that the phase coherence will 
be more important in a more specific framework which as¬ 
sumes global spiral shock waves extending across a range of 
radii in the disc. Phenomenologically, the latter model can 
be considered as a subset of the bright-spot scenario where 
the azimuthal position of individual spots coincides with a 
shock wave propagating through the disc material. 

Different assumptions stand behind the two approaches 
discussed in this work, the simplified analytic approach of 
Abramowicz et al. (1991) and our present numerical treat¬ 
ment. Apparent discrepancies between the corresponding 
values of the power-spectrum indexes, a vs. a, call for deeper 
investigation. 

Our current contribution employed energy-integrated 
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Figure 4. No phase coherence is seen in this fi^re. This is also 
indicated by formal tests: a, and p (cp. Tab. U). The resulting 
signal corresponds to a large number of spots {Ns = 200) dis¬ 
tributed over a rather broad range of radii {Q Rg < Rs < 25 Rg) 
plus the 5% noise component. Moderate inclination {Oq = 20°) 
means that variability due to orbital motion is much reduced and 
relativistic effects unimportant. We conclude that for these rea¬ 
sons the phase coherence has been suppressed. 

light curves. It should be noted that light curves resolved 
in time and energy contain more information: Standard ap¬ 
proaches of the Doppler tomography could be then used to 
determine structure on the surface of accretion discs (Karas 
&; Kraus 1996). ASCA observations offer such data (e.g. Iwa- 
sawa et al. 1996) but the resolution has not been enough yet 
to carry out this analysis. 

The author acknowledges helpful suggestions and crit¬ 
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ticipants of the Trieste-Nordic Workshop on the Bright- 
Spot Model 1996, in particular professors M. Abramowicz, 
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ported by the grants GACR 205/94/0504 and GACR 202/ 
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